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INTRODUCTION 


In  the  application  of  the  mission  analysis  required  by  MIL-A- 
008866A  (USAF)  to  fighter  and  attack  aircraft  a  problem  arises  in 
the  selection  of  the  point  in  the  sky  (velocity,  altitude,  and 
weight)  for  the  load  factor  spectrum  for  the  combat  segment  of  the 
mission.  It  can  be  shown  that  in  many  cases  important  differences 
in  the  spectrum  can  be  obtained  from  two  "reasonable"  point  selections. 

The  problem  has  been  particularly  severe  on  some  existing 
aircraft  in  that  a  ten  percent  shift  in  the  stress  spectrum  can 
produce  a  factor  of  two  change  in  life.  Therefore,  when  it  is 
considered  that  essentially  all  of  the  fatigue  damage  for  fighter 
and  attack  aircraft  is  done  in  the  combat  segment,  this  part  of 
the  mission  deserves  special  attention. 

From  an  examination  of  available  VGH  data,  it  is  evident  that 
in  both  the  air-to-air  and  air-to-ground  operations  a  fairly  wide 
variation  in  velocity,  altitude,  and  weight  is  observed.  Therefore, 
it  would  be  surprising  if  a  single  point  in  the  sky  would  provide 
an  accurate  prediction  of  the  stress  spectrum  for  a  control  point. 

This  is  even  more  evident  for  those  aircraft  which  experience  non- 
linearities  in  the  aerodynamic  data  (i.e.,  tip  stall). 

One  possible  solution  is  to  use  multiple  points  in  the  sky 
for  this  calculation.  This  can  be  effectively  accomplished  by 
taking  the  points  in  the  sky  and  their  relative  frequency  of  occur¬ 
rence  that  are  obtained  from  that  portion  of  the  fleet  that  is 
equipped  with  multichannel  recorders  (twenty  percent  of  the  fleet, 
which  is  consistent  with  current  policy,  is  believed  to  be  an 
adequate  sample).  This  can  be  done  by  taking  the  VGH  histogram 
(the  relative  frequency  of  airspeed,  normal  load  factor,  altitude, 
and  weight)  and  dividing  by  the  total  number  of  load  occurrences  to 
obtain  the  probability  that  a  load  will  occur  in  a  given  interval 
of  airspeed,  normal  load  factor,  altitude,  and  weight.  A  stress 
level  is  selected  and  a  summation  is  made  for  each  such  probability 
where  the  corresponding  stress  at  the  midpoint  of  the  interval  of 
airspeed,  etc,  is  greater  than  the  selected  stress  level.  This 
computation  produces  the  cumulative  probability  of  exceeding  a 
stress  level.  Since  the  intervals  used  for  the  data  collection 
were  not  designed  for  this  calculation,  a  provision  is  made  to 
subdivide  the  intervals  to  improve  the  accuracy  of  the  calculation. 

This  technique  is  explained  in  Section  III.  From  the  cumulative 


probability,  the  number  of  stress  exceedances  per  hour,  the  proba¬ 
bility  density  function,  and  the  stress  spectrum  can  be  obtained. 

Having  the  functions  referred  to  above  for  a  number  of  control 
points  that  is  adequate  to  cover  the  aircraft  structure  (this 
number  may  have  to  be  obtained  by  trial  and  error),  one  may  generate 
the  full  scale  aircraft  spectrum  by  assuming  that  an  arbitrary  loading 
at  the  control  points  of  the  structure  can  be  derived  from  a  linear 
combination  of  the  loading  imposed  by  balanced  load  conditions.  If 
Np  control  points  are  used,  then  Np  balanced  load  conditions  are 
used  to  represent  the  control  point  load.  The  use  of  "representative" 
balanced  load  conditions  should  provide  a  satisfactory  interpolation 
between  control  points.  These  intermediate  points  should  be  spot 
checked  against  the  true  spectrum  to  see  if  the  control  point  coverage 
is  adequate. 

Of  course,  for  a  new  design,  the  VGH  data  does  not  exist  and 
consequently  direct  application  of  this  method  is  impossible.  In 
some  cases  it  will  be  possible  to  overcome  this  difficulty  by  taking 
existing  VGH  data  from  older  aircraft  and  by  use  of  judgement  adapt  it 
to  this  procedure.  In  any  event,  the  method  should  be  applied  when  the 
proper  data  becomes  available  so  that  by  suitable  tests  and  analysis 
the  appropriate  changes  may  be  made  in  the  aircraft  life  predictions. 

One  important  application  of  this  procedure  is  fighter/attack 
tracking.  The  unusual  technique  is  to  use  the  fleet  counting 
accelerometer  data  and  compute  the  stress  for  a  single  point  in  the  sky 
that  is  believed  to  be  representative  of  the  particular  mission  flown 
(i.e.,  air-to-air  or  air-to-ground).  In  lieu  of  this  approach,  one 
could  compute  from  the  VGH  data  the  conditional  probability  of  exceeding 
a  stress  given  the  normal  load  factor.  If  this  function  were  available, 
it  would  be  possible  to  track  to  any  desired  probability  on  even  multiple 
probability  levels  depending  on  what  results  are  desired.  This  function 
can  be  generated  from  this  program  by  setting  all  occurrences  equal  to 
zero  except  those  that  fall  in  the  desired  load  factor  interval.  The 
high  positive  and  low  negative  load  factors  may  require  an  extrapolation 
from  neighboring  load  factors  because  there  may  be  too  few  data  points 
to  adequately  describe  these  functions. 

The  program  that  is  discussed  in  this  report  is  based  on  the  load 
occurrences  in  the  VGH  histogram  being  dependent  on  indicated  airspeed, 
normal  load  factor,  altitude,  and  weight.  The  stress  function  is 
based  on  the  same  quantities.  An  immediate  alternate  that  is 
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included  is  to  use  equivalent  load  factor  instead  of  load  factor. 
This  removes  the  weight  dependency  and  considerably  reduces  the 
magnitude  of  the  input.  This  option  is  included  in  the  computer 
program  described  in  the  text.  Other  alternates  that  could  be 
obtained  by  a  simple  modification  of  the  program  are  listed  as 
follows: 


VGH  data 
based  on 

1 .  Indicated  airspeed, 
normal  load  factor, 
altitude,  and  weight. 

2.  Mach  no.,  normal 
load  factor, 
altitude,  and  weight. 

3.  Equivalent  airspeed, 
normal  load  factor, 
and  weight 


Stress  function 
based  on 

Mach  no.,  normal  load 
factor,  altitude,  and 
weight. 

Mach  no . ,  normal  load 
factor,  altitude,  and 
weight. 

Equivalent  airspeed,  normal 
load  factor,  and  weight. 


The  extension  of  this  program  to  include  other  degrees  of  freedom 
for  the  aircraft  is  immediately  evident.  The  major  difficulty  is 
the  management  of  the  input  data  required  for  the  load  occurrences 
and  the  stress  function. 
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SECTION  II 


ANALYTICAL  DERIVATION  OF  THE  SPECTRUM 

The  first  step  in  the  derivation  of  the  fatigue  spectrum  is 
to  solve  for  the  stress  probability  distribution  function.  This 
requires  that  the  histogram  of  occurrences  i  n  interval  s  of  indi¬ 
cated  airspeed,  load  factor,  altitude, and  weight  be  defined.  To 
do  this  suppose  that  each  of  N  . ,  N  ,  N, ,  and  N  is  a  positive 
integer  and  vi  nz  n  w 

(1)  v^.  is  a  simple  graph  such  that  the  x-projection  of  v.. 

is  the  set  of  integers  in  [1,  Nyi  +  1]  and  if  i  is  an 

integer  in  [1,  N  .  +  1]  and  i  +  1  is  in  [1,  N  .  +  1] 

then  the  indicated  airspeed  v..(i)  is  less  than  the  in¬ 
dicated  airspeed  v^.  (i  +  1) 

(2)  n  is  a  simple  graph  such  that  the  x-projection  of  n 

zi  zi 

is  the  set  of  integers  in  [1 ,  Np  +  1]  and  if  j  is  an 

integer  in  [1,  N  +  1]  and  j  +  1  is  in  [1 ,  N  +1] 

nz  nz 

then  the  normal  load  factor  nz  (j),  is  less  than  the 

normal  load  factor  nz  (j  +  1) 

(3)  h.j  is  a  simple  graph  such  that  the  x-projection  of  h.. 
is  the  set  of  integers  in  [1,  Nh  +1]  and  if  k  is  an  in¬ 
teger  in  [1,  N^  +  1]  and  k  +  1  is  in  [1,  Nh  +  1]  then 
the  altitude  h^k),  is  less  than  the  altitude  h^(k  +  1) 

(4)  is  a  simple  graph  such  that  the  x-projection  of  w.  is 
the  set  of  integers  in  [1,  N  +  1]  and  if  m  is  an  integer 
in  [1,  N  +  1]  and  m  +  1  is  in  [1,  N  +  1]  the  weight 
w^m),  is  less  than  the  weight  w..(m  +  1) 
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Further,  suppose  Hj  is  a  simple  surface  such  that 

[v i i ( i ) »  r>z  (j),  h ^ ( k ) ,  w^m),  Hj (v^ i ( i ) ,  nz  (j),  h.(k), 
w.(m))]  is  a  point  of  H,  only  if 

1  U 

(1)  i  is  in  [1,  Nvi ] ,  j  is  in  [1,  Nn  ],  k  is  in  [1,  Nh], 

m  is  in  [1 ,  N  ]  and 
w 

(2)  H,(v..(i),  n  (j),  h.(k),  w.(m))  is  the  number  of  "load 

J  n  z.j  i  i 

occurrences"  in  the  rectangular  interval 
[vii(i),  v  1  i  ( i  +  1);  nz  (j),  nz  (j  +  1);  h^k), 

h^k  +  1);  w.(m),  w^m  +  1)] 

and  these  load  occurrences  are  assumed  to  be  uniformly 
distributed  within  the  rectangular  interval. 

The  surface  H,  is  called  the  VGH  histogram  for  v . • ,  n,  , 

0  M  z  ^ 

h^  and  wi . 

The  total  number  of  load  occurrences  included  in  the  VGH 
histogram  Hj  is 

N  N  M  N 

vi  "2  Nh  Nw 

Nt  =  X  X  X  X  Hj(v.,(i),  n  (j),  h . (k)k  w.(m)) 

z  i=l  j=l  k=l  m=l  d  11  zi  1  1 

Therefore,  by  definition,  the  probability  that  the  indicated 
airspeed,  normal  load  factor,  altitude,  and  weight  is  in  the 
rectangular  interval  [v..(i),  v..(i  +  1);  n  (j),  n  (j  +  1); 

11  1  zi  Ai 

hi  (k),  h.(k  +  1);  w^m),  w1  (nr,  +  1)]  is 


Now  suppose  that  if  i  is  in  [1,  Ny  -1]  then  the  interval 
Cv -  - ( i ) ,  Vij(i  +  1)]  is  covered  by  Nu  equal  intervals,  if  j  is  in 

li  II  V  ^ 

[1,  N  -  1]  then  [n  (j),  n  (j  +  1)]  is  covered  by  NR  equal 
nz  zi  zi  nz 

intervals,  if  k  is  in  [1,  Nh  -1]  then  [h^(k),  h^(k  +  1)]  is  cov¬ 
ered  by  NR  equal  intervals,  and  if  m  is  in  [1,  N  -1]  then  [w.(m), 
n  r  w  1 

w.  (m  +  1)]  is  covered  by  N  equal  intervals. 

1  W 

Since  it  was  supposed  that  the  load  occurrences  are  within  the 
rectangular  interval  [v.^i),  v..(i  +  1);  nz  (j), 

nz  (j  +  1);  h.(k),  h^(k  +  1);  w. {m),  w^(m  +  1)]  then  the  probability 

that  the  indicated  airspeed,  normal  load  factor,  altitude,  and 
weight  is  in  the  rectangular  interval 


is 


v. i (i  +  1)  -  v . . ( i ) 
[v4<<1).  v„0)  +  -U-  " 


n 


n 


n  (j),  n  (j)  + 
i  zi 


N 

Vi 

n  (j  +  1)  -  n  (j) 
zi  zi 


N' 


h.(k  +  1) 

h,(k),  h.(k)  +  -J - 


z 

hjOO 


N 


w. (m  +  1)  -  w.(m) 
wi(m),  w.(m)  + - - - - ] 


N, 


w 


a  Hj(v.(i),  nz(j),  h(k),  w(m)) 

P-jO»  J»  k,  m)  =  ■■  b  d  d  b 

w  at  at”  at*' 


(1)  Now  suppose  that  is  a  simple  surface  such  that  the 
x-y  projection  of  V..  is  the  set  of  integers  in  the  rectangular 
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D 

interval  [1 ,  N  +  1;  1,  N  ]  and  if  i  and  i  +  1  are  ingers 
vi  vi 

n 

in  [1,  Nv .  +  1]  and  iR  is  an  integer  in  [1,  ]  then 

V-  V  >  <*,i«  ♦’>  - 

vi 

(2)  Nz  is  a  simple  surface  such  that  the  x-y  projection  of 

N  is  the  set  of  integers  in  the  rectangular  interval 
Z  R 

[1 ,  Nn  +  1 ;  1 »  Nn  ]  and  if  j,  and  j+  1  are  integers  in 

n 

[1 ,  +1]  and  jr  is  an  integer  in  [1,  Nn  ]  then 

jp  -  0.5 

Nz(j,  JR)  •  nz  (j)  +  (-^ - )  (n  (j  +  1)  -  n  (j)) 

i  N  i  i 

n 

z 

(3)  H  is  a  simple  surface  such  that  the  x,  y  projection  of  H  is 

the  set  of  integers  in  the  rectangular  interval  [1 ,  N.  +  1 ; 

R  n 

1,  Nh]  and  if  k  and  k  +  1  are  integers  in  [1,  Nh  +  1]  and 

kR  is  an  integer  in  [1,  N*-,  then 

k  -  0.5 

H(k,  kR)  =  h.(k)  +  (-^ - )  ( h .  ( k  +  1)  -  h.(k)) 

1  % 

(4)  W  is  a  simple  surface  such  that  the  x,  y  projection  of  W 
is  the  set  of  integers  in  the  rectangular  interval  [1, 

D 

Nw  +  1 ;  1,  Nw]  and  if  m  and  m  +  1  are  integers  in  [1,  Nw  +  1] 

D 

and  mp  is  an  integer  in  [1,  N  ]  then 

m  -  0.5  w 

W(m,  mR)  =  w.  (m)  +  (~^ - )  (w.  (m  +  1)  -  (m)) 

Nw 


The  assumption  is  made  that  the  stress  at  a  point  in  the  structure 
depends  only  on  the  indicated  airspeed,  normal  load  factor,  altitude 
and  weight.  Therefore,  if  it  is  supposed  that  each  of  a  and  N  is  a 
positive  integer  such  that  a  is  in  [1,  N  ]  p 
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and  ¥a  is  a  simple  surface  such  that  (V^ ( i , iR) ,  Nz(j,jR),  H(k,kR), 

W(m,mR),  4'a(Vi ( i »iR) »  Nz(j,jR),  H(k,kR),  W(m,mR))  is  a  point  of 

y3  only  if  i  is  in  [1 ,N ,  +  1],  iD  is  in  [1,N^  ],  is  in 

vi  K  vi 

[1»NW  +  1],  mR  is  in  [1,N*]  and  H-a(Vi(i,iR),  Nz(j,jR),  H(k,kR), 
W(m,mR))  is  the  stress  for  the  ath  control  point  corresponding  to 
the  indicated  airspeed  V^(i,iR),  the  normal  load  factor  Nz(j,jR), 
the  altitude  H(k,kR),  and  the  weight  W(m,mR). 

a  A 

The  surfaces  v  and  Pj  are  used  in  the  calculation  of  the 

cumulative  probability  of  exceeding  a  given  stress  as  follows: 
Suppose  that  N  is  a  positive  integer  and  is  a  uniformly  in- 

a 

creasing  sequence  with  x-projection  [1,N  ]  and  $  is  a  simple 

rL 

surface  such  that 

(1)  $  (i ,j ,k,m,iR,  j R , kR ,mR )  -  Pj(i»j»k,m) 

if  ¥a(Vi(i,iR),  Nz(j,jR),  H(k,kR) ,  W(m,mR))  >  rL(b) 

(2)  $  ( i » j » k,m, iR, jR , kR ,mR )  -  0 

if  ¥a(Vi(i,iR),N2(j,jR),H(k,kR),W(m,mR))  $  rL(b) 


rL(b)  is 


Therefore,  the  probability  that  the  stress  is  greater  than 

R  nR 


N, 


Vi 


N 


Pya(  rL(b) )  =  X  X 
i=l  j=l 

$  ( i » j t k ,m, iR , jR , kR ,mR ) 


Nh  Nw 

X  X 

k=l  m=l 


N 


Vi 


rr  ,,r 

"z  "h 

X  X 


N 

w 

X 


iR=1  jR  1  kR=l  mR-l 


The  probability  density  function  PQ  is  the  derivative  of  the 

cumulative  probability  function  P^a.  This  derivative  is  computed 
as  follows:  Suppose  a  is  an  integer  in  [1,N  ]  and  that  ^a  is  a 
simple  graph  with  x-projection  the  interval  [l,Nr]  such  that 
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(1)  if  b  is  an  integer  in  [1,  Nr  ]  then  ca(b)  =  P  a (b )  and 

L 

(2)  if  c  is  a  number  in  [b,  b  +  2]  there  exists  a  u-j ,  u^, 
and  such  that  ?a{c)  +  u.|C2  +  U2C  +  where  ,  u^» 
U3  are  determined  from  the  equations 


Ca(b) 

b2 

b 

1 

U1 

U2 

ta(b  +  1) 

(b  + 

l)2  (b 

+  1) 

1 

C3(b  +  2) 

(b  + 

2  )2  (b 

+  2) 

1 

u3 

Therefore 

(1 )  if  b  =  1 

V1>  =  2“> 

rL0) 

+  U2 

V2) =  2  “> 

rL(l ) 

+  u2 

(2)  if  b  is  in  [2,  Nr  -  3] 

L 

Vb  +  ”  =  2  U1  rLlb  +  ’>  +  u2 


(3|  if  b  »  Nr  -2 

V(V,)  =  2“'  rL(\‘1)+“2 
VV  =  2  u>  rL(V  *  “2 

The  next  step  in  the  derivation  of  the  fatigue  loading  spectrum 
is  to  determine  the  stress  and  the  frequency  of  that  stress  in  the 
spectrum.  This  is  done  by  an  indirect  process  as  shown  below. 

Suppose  that  the  fatigue  test  spectrum  is  to  be  composed  of  N  cycles 
at  M  stress  levels.  Further,  suppose  that  a  is  a  positive  integer  in 
[1,  Np]  and  Sa  is  a  sequence  of  M  numbers  such  that  sa(i)  and 

sa(j)  are  members  of  Sa  only  if  0  <  sa(i)  <  sa(j)  <  1  and  i  <  j. 


Therefore,  each  member  of  Sa  corresponds  to  an  ordinate  of  the 
graph  1 -Pya .  The  M  abscissas  corresponding  to  these  M  ordinates 

are  defined  as  the  M  stress  levels  of  the  spectrum  for  the  ath 
control  point.  The  graph  1-P  a  is  known  at  N  points.  Conse- 

rL 

quently,  an  approximation  to  1-Pa  must  be  found  in  order  to  com¬ 


pute  the  spectrum  stress  levels.  Suppose  a  is  a  simple  graph  with 
x-projection  the  interval  [r,  (1),  r,  (N ^  )]  and  if  k  is  in  [1  ,N  ] 

L  L  rL  rL 

then  8  (rL(k))  =  1 -P^a ( rL ( k ) ) .  Further,  suppose  that  if  i-1,  i, 

and  i  +  1  are  in  [1,  ],  6L  is  rL(k  +  1)  *  r^(k),  and  x  is  in 


[-<$L,  <5l]  then 

e(x  ♦  rL(k))*  D  (f-)2] 


r  o  i 


o 

i 

7 


'&(  rL(k-l  ))N 

I  s(rL(k)) 

i 

|  e(rL(k+l))  j 


\  >  \  > 

It  follows  then  that  if  i  is  in  [1,M)  there  exists  an  inte¬ 
ger  k  such  that  pya(rL(k-l))  <s  (i)  <  P^a  (rL(k))  and  a  number  x 

such  that  8 (x  +  rL { k ) )  =  sa ( 1 ) .  The  number  x  is  obtained  from  a 
solution  of  a  quadratic  equation  and  x  +  rL(k)  is  the  stress  cor¬ 
responding  to  sd(i). 


The  fraction  of  the  N  cycles,  n^,  that  are  associated  with 
the  ith  stress  level  is  defined  as  follows: 

n,  . 

.  imi  -  s(<-i)  , , ,  <  „ 

^ . , .  im  *  stf-ij) 

It  follows  then  that  if  i  is  in  [1,M]  and  if  the  sequence  Sais 
used  for  each  of  the  control  points  then  there  will  be  an  equal 
number  of  loading  cycles  for  the  ith  load  level  for  each  of  the 
control  points. 


The  final  step  Is  to  determine  a  set  of  coefficients  which 
when  multiplied  by  the  stresses  corresponding  to  balanced  load 
conditions  for  the  aircraft  will  produce  the  desired  stress  levels 
at  the  aircraft  control  points.  Suppose  a  is  a  positive  integer 
In  [1,M],  b  is  a  positive  Integer  In  [1,  N  ],  and  c  is  a  positive 

Integer  in  [1,N  ].  Therefore,  'f  A*b  is  the  stress  for  the  ath 

load  level  at  the  bth  point  In  the  sky  anci  the  cth  control  peine 
and  r*  Is  the  stress  desires  In  the  fatigue  test  for  che  ath  load 

level  and  the  ctn  control  point  then  there  exists  a  set  of  coet- 
flcents  a*b  such  that  i|  *  aab. 
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SECTION  III 


DESCRIPTION  OF  THE  COMPUTER  PROGRAM 

1  NOTATION 

The  right  hand  side  of  the  following  relations  are  defined  in 
Section  II. 

NT 4 21  =  Ny 

NT422  =  N 

n2 

NT423  =  Nh 
NT424 * N 

w 

NT  •  Nt 

PJT  =  Pj 

NRVI  =  NR 
V1 

NRNZ  =  NR 
nz 

NRH  =  NR 

NRW  =  NR 
w 

vn  - 

NZI  =nz 

HI  =  b. 

WI  =  w . 

VI  =  vi 

NZ  =  Nz 

H  =  H 
u  =  u 

NPS  =  Np 

PPSI  =  P^a  (The  a  is  not  explicitly  identified  in  the  program) 
PDPSI  =  PQ  (The  a  is  not  explicitly  identified  in  the  program) 

PS  =  A 
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PLD  =  r 

ALPHA  =  a 

PV I  =  v . .  ( N  +1) 
u  v. 

FNZ  =  n  (N  +  1 ) 
zi  nz 

FH  =  h.(Nh  +  1) 

FW  *  w.(Nw  +  1) 

FACT0R  -  Stress  scaling  factor.  FACT0R  =  1  unless  otherwise 
specified. 

H0URS  -  The  number  of  hours  of  data  in  the  VGH  histogram 
PSIL  =  rL 

ARLAN  =  sa 
DELTA  =  6l 

APDPSI ( I )  =  1.0  -  PPSI(I) 

PSILL  (I)  -  The  stress  level  that  is  the  abscissa  of  the 
point  of  PPSI  whose  ordinate  is  AREAN(I) 

FRAC(I)  -  The  fraction  of  the  total  number  of  cycles  in  the 
spectrum  that  correspond  to  PSILL 

NPSIL  =  Nr 

L 

N PSILL  =  M 

EXCEED ( I )  -  The  number  of  exceedances  per  hour  of  the  stress 
PSIL(I) 

NZERO  -  Control  number  to  zero  the  input  numbers  at  the  start 
of  a  run  and  then  prevent  them  from  being  zeroed 
between  cases 

NPSCT  -  Control  number  for  counting  the  number  of  control 
points  for  which  a  spectrum  has  been  computed  in  a 
single  run 

2  INTERPOLATION  PROCEDURE 

Since  the  stress  is  initially  calculated  for  only  a  finite 
set  of  points  on  the  stress  surface,  an  assumption  must  be  made 
to  determine  the  stress  for  a  given  indicated  airspeed,  normal 
load  factor,  altitude,  and  weight.  Specifically,  the  problem  may 
be  expressed  as  follows:  Given  that  NT421 ,  NT422,  NT423,  NT424, 

NRVI,  NRNZ,  NRH,  and  NRW  is  a  positive  integer  and  I  is  in  [1.NT421], 
J  is  in  [1.NT422],  K  is  in  [1.NT423],  M  is  in  [1,NT424],  IR  is  in 
[l.NRVI],  JR  is  in  [1,NRNZ],  KR  is  in  [l.NRH],  MR  is  in  [1 ,NRW]  and 
a  is  in  [l,Np]  it  is  required  to  create  an  approximation  in  the  form 


1 

i 

A 

i 
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£*(I ,J,K,M,IR,JR,KR,MR)  = 

=a(VI(I,IR),NZ(J,JR),H(K,KR),W(M,MR)) 
for  the  stress  as  expressed  by 

,j,a(I ,  J,K,M,IR,  JR,KR,MR)  = 

fa(VI{I,IR),NZ(J,JR),H(K,KR),W(M,MR)) 

a  4  3 

where  =  is  a  ruled  surface  based  on  2  =16  points  of  f  .  The 

method  of  choosing  these  16  points  and  the  calculation  of  the 
stress  approximation  is  described  below. 

The  first  step  is  to  define  the  function  TABLE  which  contains 
the  projections  and  ordinates  of  the  fa  surface. 

Suppose  each  of  NTAB1 ,  NTAB2,  NTAB3,  and  NTAB4  is  a  positive 
integer  and  that 

NN1 2  =  NTAB1  +  NTAB2 

NN1 3  =  NN12  +  NTAB3 

NN1 4  =  NN13  +  NTAB4 

NP  =  NTAB1  •  NTAB2  *  NTAB3  *  NTAB4 

NF  =  NN1 4  +  NP 

Further,  suppose  that  TABLE  is  a  simple  graph  such  that  the 
x-projection  of  TABLE  is  the  set  of  integers  in  the  interval  [l.NF] 
and  each  of  II,  12,  13,  and  14  is  a  positive  integer. 

Also, 

(1)  if  II  and  II  +  1  are  in  [1.NTAB1]  then  the  indicated 
airspeed  TABLE  (II)  is  less  than  the  indicated  airspeed 
TABLE  (II  +1) 

(2)  if  12  and  12+1  are  in  [NTAB1  +  1,  NN12]  then  the  normal 
load  factor  TABLE  (12)  is  less  than  the  normal  load  factor 
TABLE  (12+1) 

(3)  if  13  and  13  +  1  are  in  [NN12  +  1,  NN13]  then  the  altitude 

TABLE  (13),  is  less  than  the  altitude  TABLE  (13  +  1) 

(4)  if  14  and  14  +  1  are  in  [NN13  +  1  ,  NN14]  then  the  weight 

TABLE  (14)  is  less  than  the  weight  TABLE  (14+1) 

(5)  if  II  is  in  [1 ,NTAB1 ] ,  12  is  in  [NTAB1  +  1,  NN12], 

13  is  in  [NN12  +  1,  NN13],  14  is  in  [NN13  +  1,  NN14] 
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and  n  is  in  [NN14  +  1,  NF]  and  is  equal  to  NN14  + 

(14  -  NN13  -  1)  *  NTAB3  •  NTAB2  *  NTAB1  +  (13  -  NN12  -  1) 
-  NTAB2  *  NTAB1  +  (12  -  NTAB1  -  1)  *  NTAB1  +  II  then  the 
stress  TABLE  (n)  is  the  stress  that  corresponds  to  the 
indicated  airspeed  TABLE  (II),  the  normal  load  factor 
TABLE  (12),  the  altitude  TABLE  (13)  and  the  weight  TABLE 
(14). 

The  positive  integers  II,  12,  13,  and  14  are  determined  as 
follows:  A  search  is  made  for  the  integer  i  that  will  determine 
the  smallest  number  TABLE  (i)  that  equals  or  exceeds  VI(I,IR).  If 
i  =  1  satisfies  this  requirement  then  II  is  set  equal  to  2.  If  i 
is  in  [2,  NTAB1]  then  II  is  set  equal  to  i .  If  no  i  can  be  found 
in  [2,  NTAB1]  then  II  is  set  equal  to  NTAB1 .  A  search  is  made  for 
the  integer  j  that  will  determine  the  smallest  number  TABLE  (j) 
that  equals  or  exceeds  (NZ(J,JR).  If  j  =  NTAB1  +  1  then  12  is  set 
equal  to  NTAB1  +2.  If  j  is  in  [NTAB1  +  2,  NN12]  then  12  is  set 
equal  to  j.  If  no  j  can  be  found  to  satisfy  the  requirement  then 
12  is  set  equal  to  NN12.  Also,  a  search  is  made  for  the  integer  k 
that  will  determine  the  smallest  number  TABLE  (k)  that  equals  or 
exceeds  H(K,KR) .  If  k  =  NN12  +  1  then  13  is  set  equal  to  NN12  +  2. 
If  k  is  in  [NN12  +  2,  NN13]  then  13  is  set  equal  to  k.  If  no  k 
can  be  found  in  [NN12  +  2,  NN13]  then  13  is  set  equal  to  NN13.  A 
final  search  is  made  for  the  integer  m  that  will  determine  the 
smallest  number  TABLE  (m)  that  equals  or  exceeds  W(M,MR).  If  m  = 
NN13  +  1  then  14  is  set  equal  to  NN13  +2.  If  m  is  in  [NN13  +  2, 

NN14]  then  14  is  set  equal  to  m.  If  no  m  can  be  in  [NN13  +  2, 

NN14]  then  14  is  set  equal  to  NN14. 

The  next  step  is  to  identify  the  integers  required  for  the 
final  calculations. 

With 

NP12  =  NTAB1  •  NTAB2 

NP13  =  NP12  '  NTAB3 

these  are: 

N2222=  NN14  +(I4  -  NN13  -  1)  *  NP13  +  (13  -  NN12  -  1)  * 

NP1 2  +  (12  -  NTAB1  -  1)  *  NTAB1  +  II 

N1222  =  N22  -  1 

N2122  “  N222  -  NTAB1 

Nil 22  =  N2122  -  1 

N2212  =  N2222  -  NP12 
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.1 

i 


N1212 
N2112 
N1112 
N2221 
NT  221 
N2121 
Nil  21 
N2211 
N1211 
N2111 
Nil  1 1 


N2212  -  1 
N2212  -  NTAB1 
N2112  -  1 
N222  -  NP13 


N2221 

N2221 

N2121 

N2221 

N2211 

N2211 

N2111 


-  1 


NTAB1 

1 

NP1 2 
1 

NTAB1 

1 


Therefore,  if 


XI  RAT  = 
X2RAT  = 
X3RAT  = 

X4RAT  = 

then 

AMP111 

AMP211 

AMP121 

AMP221 

AMP112 

AMP 21 2 

AMP122 

AMP222 

AMP11  = 
AMP12  = 
AMP22  = 


« 


rm 

NZ(J,JR)  - 

TABLE  (I 


TABLE  (11-1) 

-■  md (tpit 

TABLE  (12-1) 


(12)  -  TABLE  (12-1 


TABLE  (13-1) 
-  TABLE 


[FIT 


H(K,KR)  - 
TABLE  (13) 

W(M,MR)  -  TABLE  (14-1 

TABirntr  "  "" 


-■We  (TC-TT 
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TABLE (N1 1 11 ) 
TABLE(N1 21 1 ) 
TABLE ( N1 1 21 ) 
TABLE(N1 221 ) 
TABLE  Nil  12 
TARLE(N1212 
TABLE (Nil  22 ) 
TABLE(N1 222) 


XI  RAT (TABLE(N211 1 ) 
X1RAT (TABLE(N221 1 ) 
XI  RAT (TABLE (N21 21 ) 
XI  RAT (TABLE (N2221  ) 
X 1  RAT (TABLE (N21 1 2 ) 
X1RAT  TARLE  N2212) 
X1RAT(TABLE(N2122) 
X1RAT(TABLE(N2222) 


TABLE ( N1 1 1 1 ) ) 
TABLE (N1 21 1 )) 
TABLE ( N1 1 21 )) 
TABLE(N1 221 ) ) 
TABLE  Nil 12)) 
TARLE(N1212) 
TABLE (Nil  22)) 
TABLE (N1 222) ) , 


AMP111 
AMP1 1 2 
AMP122 


AMP1  =  AMP11  + 
AMP 2  =  AMP12  + 


+  X2RAT (AMP21 1  -  AMP 111 ) 

+  X2RAT (AMP21 2  -  AMP112) 

+  X2RAT (AMP222  -  AMP122), 

X3RAT (AMP21  -  AMP 11 ) 

X3RAT (AMP22  -  AMP! 2), 


Z  (K.J.K.M.IR.JR.KR.MR)  =  (AMP!  +  X4RAT(AMP2  -  AMP1 ) ) • FACT0R 

It  is  seen  that  the  sixteen  points  on  the  y3  surface  are  reduced 
to  eight  points  on  the  =  surface  by  an  interpolation  on  the  indi¬ 
cated  airspeed.  The  eight  points  are  reduced  to  four  points  on  the  = 
surface  by  an  interpolation  on  the  normal  load  factor.  Next,  the 


four  points  are  reduced  to  two  points  on  the  =d  surface  by  an  in¬ 
terpolation  on  the  altitude,  and  finally  these  two  points  are  re¬ 
duced  to  the  desired  stress  by  an  interpolation  on  the  weight. 

Note  that  the  number  FACT0R  is  used  to  scale  the  calculation 
made  in  the  table  look  up  routine. 

3  COMPUTER  FLOW  DIAGRAM  AND  PROGRAM 

The  computer  routine  was  coded  in  FORTRAN  Extended  Language 
with  the  main  program  and  subroutines  arranged  as  follows: 


SPECA  Program  (See  Appenoix  A) 


- *:  TABT  i 


iPAGEHDk 


MAIN 


GUIDE  k- 


INPUT 


INAID 


CALC 


f= 


-LDLVL 


PRINTRi 


jLDC0EFf*- 


TABQ 


«*GAUSS2 


-GAUS21 


MAIN  -  Main  Program  -  Sets  NZER0  and  NPSCT  to  zero  and  trans¬ 
fers  program  control  to  GUIDE 
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GUIDE  -  Subroutine  -  Initially  zeros  input  and  output  numbers 
and  after  first  case  zeros  output  numbers  before  the 
calculations  are  performed.  GUIDE,  tne  main  control¬ 
ling  subroutine,  transfers  control  to  INPUT,  CALC, 
lDlVL,  PRINTR,  and  LDC0EF  in  turn. 

INPUT  -  Subroutine  -  Reads  in  all  input  data  including  the 
VGH  histogram  and  the  stress  table.  There  are  two 
formats  for  reading  in  floating  point  numbers  and 
three  formats  for  reading  in  fixed  point  numbers. 

Tne  details  of  the  data  input  are  discussed  later  in 
this  section. 

INAID  -  Subroutine  -  Called  by  INPUT  and  has  the  purpose  of 
writing  out  certain  input  data. 

(1)  NRV I ,  NRNZ ,  NRH,  NRw 

(2  FACT0R 

(3)  PSIL 

(4)  AREAN 

(5)  PS 

(b)  Stress  table 

(7)  V3H  histogram  table 

(8)  FVI ,  FNZ,  FH,  FW 

Also,  INAID  sets  NZER0=1  for  control  of  data  handling 
in  GUIDE 

PAGCHD  -  Subroutine  -  Writes  out  page  heading  inducing  run 
identification,  date  and  page  number 

CALC  -  Ssoroutme  -  Computes  PPSI  and  PDPSI 

LDLVL  -  Subroutine  -  Computes  PSIL  and  F.IAC 

LDC0EF  -  Subroutine  -  Computer  ALPHA 

TABQ  -  Subroutine  -  Called  from  CALC  to  perform  the  interpo¬ 
lation  discussed  in  Section  IV,  13,  that  computes  the 
stress  corresponding  to  a  given  indicated  airspeed, 
normal  load  factor,  aioicude,  and  weight. 

TABT  -  Subroutine  -  called  from  CALC  as  an  alternate  to  TABQ 
for  the  interpolation  to  compute  the  stress  correspon¬ 
ding  to  a  given  indicated  airspeed,  equivalent  normal 
load  factor,  and  altitude. 
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GAUSS2  -  Subroutine  -  Called  from  CALC  to  solve  the  simu.tun- 
eous  equations  that  are  required  to  pass  second  order 
equations  through  the  points  of  PI-SI  so  that  the  dif¬ 
ferentiation  vor  PDPS1  can  be  perforned.  The  sub¬ 
routine  uses  the  Gauss-Jornan  method  for  solving  the 
sets  ov  simultaneous  equations. 

GAUS2I  -  Subroutine  -  Called  from  LUC0EF  and  is  used  to  solve 
the  set  of  equations  PLD(l)  =  P3(K.,J)  '  ALPHA (J ) . 

Ihis  subroutine  is  identical  to  GAU3S2  except  for  a 
U1MENSI0N  statement  change. 

PkiNTR  -  Subroutine  -  Called  from  GUIDE  to  write  out  computed 
output  data.  In  particular,  PR1NTK  prints 

(1)  PPSI 

(2)  tXCEEU 

13)  H0URS 

(4)  PDPSI 

(b)  PSILL,  FRAC 

4  EQUIVALENCE  TABLES 

The  technique  that  has  been  used  in  coding  this  routine  is  to 
place  all  input  and  output  numbers  In  blank  common.  All  input  and 
output  floating  point  numbers  are  called  parameters  and  are  con¬ 
tained  in  P  (dimensioned  10,U00).  All  input  ano  output  fixea  point 
numbers  are  called  integers  and  are  contained  in  NTEGER  (dimensioned 
I0U).  To  make  the  program  more  easily  interpreted,  EQUIVALENCE 
statements  a'-e  used  tc  provide  the  ?  and  NTEGER  numbers  w'th  more 
recognizable  names.  The  5PECA  program  parameter  and  integer  sables 
are  given  below. 

PARAMETER  EQUIVALENCE  TABLE _ 


I  P  1  Dimension 

Term 

1  P  -- 

Dimension 

Term 

T - - 

nr 

FMN.fVI 

1  1  2  J 1 

(100) 

APDPSI ( 1 ; 

2 

(i) 

FNZ 

:  I3cn 

1 

APDPSI (100, 

3 

(i) 

FH 

*  1301 

(IOC) 

PSIL(I) 

4 

a) 

FW 

•  1400 

PSIL(IOO) 

5 

a) 

FACT0R 

1  MOI 

1  (100) 

FRAC ( 1 ; 

6 

m 

H0URS 

1  1S00 

1 

FKAC(I0U:  1 

1  SOI 

:  (100,25) 

PLDS( 1,1)  1 

• 

,  4000 

i 

PLDS(i 00,25)  ; 

ALPHA ( 1 )  i 

• 

i  4001 

!  U5) 

• 

4025 

i 

ALPHA (25)  ! 

i 
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5  INPUT  DATA 


AH  of  the  input  data  described  below  is  read  into  the  pro¬ 
gram  by  means  of  the  subroutine  INPUT,  INPUT  is  a  general  purpose 
subroutine  for  reading  data  from  cards.  For  this  program,  the  full 
capabilities  of  INPUT  are  not  required  and  consequently  there  will 
be  some  zeros  in  the  input  that  serve  to  bypass  certain  options. 

The  following  deck  arrangement  is  recommended: 


1415  Format 


IDENnNPFl 


0  ,  0 
_ I _ 


0  NTI4.HTW4: MONTH 


DAY 

YEAR 

NPSIL 

npsillI 

NPS  |  21 

715  Format 


T 


INRVI  NRNZ  NRH 


NKW 


NIB 


_ 72H  Format 

Run  Description 


72H  Format 

i 

Run  Description 


315 


6E10.3 

i  Format 

FVI  1 

I 

FNZ  | 

FH 

1  FW 

'  — . .  - ,  ,i 

FACT0R  ] 

|  H0UKS 

___3I5  Format 
jlOl 


100  + 
NPSIL 


m 
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6E10.3  Format 


PSIL(l)  -  PSIL{NPSIL) 


315  Format 
20llZ00  +  |  , 
^Ul | NPS1LL  [  1 


6E10.3  Format 


AREAN(l)  -  AKEAN(NPSILL) 


If  NPS  >  1  go  to  (a) 
If  NPS  =  0  go  to  lb) 


(a) 


301 


315  Format 

T00+ 


NPS 


I 


6E10.3  Format 


PS( 1,1)  -  t PS ( NPS , 1 ) 


315  Format 


326 


325  + 
NPS 


1 


6E10.3  Format 


PS(1,2)  -  PS(NPS,2) 

I 

i - - -  - - -  ■  -  -  ■  -  - 
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315  Format 


351 


.350  + 
:NPS 


1 


6E10.3  Format 


PS(1 ,3)  -  PS(NPS,3) 


315  Format 


300  + 

- 1 

,25(NPS 

(NPS-I) 

1 

ku _ 

•25+NPS 

6EI0.3  Format 


PS ( 1 ,NPS)  -  PStNPS.NPS ) 


(b)  If  NTH  >  0  go  to  (c)  to  read  N 1 14  tablets) .  For  the  first 
run  in  a  computer  input  the  stress  table  and  the  VGH  histo¬ 
gram  table  must  be  read.  Subsequent  runs  may  require  no  new 
tables  (NTI4  =  0),  one  new  table  (NTI4  =  I),  or  two  new  tables 
t  NT  1 4  =  2). 


If  NTH  =  0  go  to  1 9 ) 


5110  1 

rormat 

rn 

NTABl 1 

NTAB21 

NTAB3 

I  NTAB4 

(Stress  table  control  cards) 
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6E1Q.3  Format 


TABLE(1 )  -  TABLt(NTABI) 

(indicated  airspeeds  for  the  stress  table) 


_ 6E10.3  Format _ 

TABLE (NTAB1+1 )  -  TABLE (NN1 2 ) 

(normal  load  factors  for  the  stress  table) 


_ _ 6E1Q.3  Format _ 

TABLE (NNJ 2+ 1 )  -  TABLE (NN 13) 
(altitudes  for  the  stress  table) 


If  NTB  =  I  go  to  (d) 
If  NTB  =  2  go  to  (e) 


(d ) _ _ 6E10.3  Format _ 

TABLE (NN 1 3+ I )  -  TABLb(NN14) 

|  (weights  for  the  stress  table) 


_ _ _ 6EIQ.3  Format _ 

TABLE (NN1 4+ 1 )  -  TABLE (NF ) 

(stress  amplitudes  for  the  stress  table) 

(see  Section  3.2  for  ordering  of  these  entries) 


go  to  (f) 


(e)  E10.3  Format 

WTTB3  ' 
(ref,  weight)  • 


_ 6E10.3  Format _ 

1ABLE(NN13+1 )  -  TABLt(NF) 

(stress  amplitudes  for  the  stress  table) 

(see  Section  3.2  for  ordering  of  these  entries) 


(f) 


5110  Format 


NT421 

NT422 

NT423 

NT424 

(VGH  histogram  control  cards) 
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6E10.3  Format 


V 1 1  ( 1 }  -  V  1 1  (NT 4 21 ) 

(indicated  airspeeds  for  VGH  histogram  table) 


NZi  ( 1  )  -  NZI (NT422 ) 

(normal  load  factor  for  VGH  histogram  table) 


L 


HI ( 1 )  -  HI (NT423) 

(altitudes  lor  VGH  histogram  table) 


;  Will;  -  W1 (NT424 ) 

j  (weights  for  VGH  histogram  table) 


:  a ( V  1 1  ( 1  ) , NZ 1  ( 1  )  *  H 1 ( I ) ,  W  i  ( 1 ))  - 

i a  (V 1 1  (NT42 1  )  .NZI (NT422) ,Hi (NT423,W1 (NT424) ) 
(load  occurrences  in  VGh  histogram  table) 

(see  discussion  below  for  ordering  of  these  entries) 


(g)  END  0F  FILE 

The  first  card  contains  14  fixed  point  (integer) 
numbers  arranged  in  lb  fields.  1  hese  14  entries  in 
order  on  this  card  are 


(1) 

1UENT 

-  run  number 

(2) 

NPF1  = 

•* ,f  %  = 1 

= 

3  *  vf  "p  ' 

1 3) 

0 

(4) 

0 

(5) 

0 

(6) 

NT  14  - 

the  number  of  quadruple  tables  to  be  read 

this  count  the  stress  table  and  the  VGH  histogram 
table  are  each  considered  quadruple  tables.) 


2b 


(/) 

NlW4  -  1  lor  print  of  quadruple  tables 
=  0  otherwise 

(B) 

MONTH  -  month  m  date  for  page  (leading 

, 

is; 

DAY  -  day  in  date  for  page  heading 

> 

(10) 

YtAR  -  year  in  date  for  page  heading 

(11) 

NPSiL  -  n 

*  1 

(12) 

NFS  ILL  •  K 

(13) 

NPS  -  1  he  number  of  control  points  if  N  •»  1 
if  N  =  1  p 

D 

1  .  NPS  -  0 

(14) 

r 
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The  second  card  contains  seven  fixed  point  numbers  arranged 


in 

15  fields.  In  order  these  entries  are 

(i) 

NKV I 

(2) 

NRNZ 

-  nr1 

n 

(3) 

NRH  = 

(4) 

NRW  = 

nr 

W 

(5) 

0 

(6) 

0 

(7) 

NTB  = 

1  if  the  load  occurrences  in  the  VGii  histogram  de 

pend  on  indicated  airspeed,  normal  load  factor, 
altitude,  and  weight. 


NftJ  =  2  if  the  load  occurrences  in  the  vGh  histogram  de¬ 
pend  on  indicated  airspeo-a,  equivalent  normal  load 
factor,  and  altitude 

The  third  and  fourth  cards  conta  in  a72Hficld  each  for  the  pur¬ 
pose  of  run  description,  etc. 

The  fifth  card  contains  1,  6,  and  1  in  15  fields 

The  sixth  caro  contains  six  floating  point  numbers 
arranged  in  E10.3  fields.  These  six  numbers  are  placed  in 
the  following  order: 

(1)  FVI  =  v.'N,  +1) 
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(2; 

FNZ  =  n  (H  +1) 

nz 

(3) 

FH  =  h(Nh+l) 

(4) 

FW  =  w(N  +1) 

w 

(5) 

FACT0R  -  stress  scaling 

factor 

(6) 

H0URS  -  number  of  hours 

of  data  in  the  VGH  histogram 

The  seventh  card  contains  the  three  fixed  point  numbers  101, 
iOC  +  NPSiL,  1  in  order  in  15  fields.  NPSIL  must  not  exceed 
100. 

The  nexc  c<ira(s)  contain(s)  the  numbers  PSIL(l)  throuqh  PSIL(NPSIL) 
in  LI 0.3  fields,  six  numoers  per  card. 

The  next  entry  contains  the  fixed  point  numbers  201,  200  + 

NPSILL ,  1  in  order  in  15  fields.  NPSILL  must  not  exceed 
100. 

Following  this  card  the  floating  point  numbers  AkLAN(l) 
througn  AREAN(i'IP$I!.),  arranged  in  E10.3  fields,  six  num¬ 
bers  per  card,  are  entered. 

It  UPS  -  0  then  the  PS  matrix  is  omitted  from  tne  input  deck. 

if  NP5  1  then  the  PS  matrix  is  placed  next  in  the  input 
neck.  PS  is  dimensioned  (25,25;  and  is  equivalencea  to 
P  sucn  that  -  PS(1,1).  Therefore,  it  follows  that 

P(300+flPS)  PS(UPS,1),  P(32o)s  PS(1,2),  and  P(301+25(NPS-1 )) 

PS ( 1 , UPS ) .  Consequently  the  NFS  blocks  of  aata  are  read 
i n  as  ro i 1 ows : 

First  block  - 

Ine  first  card  contains  the  fixed  point  numoers  301, 

300+NPS,  1  arranged  in  15  fields. 

The  next  entries  are  the  floating  point  numbers  PS(1,1) 
through  PS(NFS,1)  in  E10.3  fields,  six  numbers  per 
card . 

Second  block  - 

The  first  card  contains  the  fixed  point  numbers  326, 

325+NPS,  1  arranged  in  15  fields. 
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The  next  entries  are  the  floating  point  numbers  PS (1,2) 
through  PS(NPS,2)  in  E10.3  fields,  six  numbers  per  card. 


Nr'Sth  block  - 

The  first  card  contains  the  fixed  point  numbers  301+ 

25 (HPS-1  ),  300+(NPS-l ) (25)  +  NPS  arranged  ,n  15  fields. 

The  next  entries  are  the  floating  point  numbers 
PS(1,NPS)  through  PS(NPS,NfS)  in  EiO.3  fields,  six 
rumbers  per  card. 

The  remaining  entries  are  the  stress  table  and  the  VGH 
histogram  table.  These  entries  are  prepared  as  follows: 

If  1  =  1  then  the  entry  is  tne  stress  table  where  there  are 
IIP  =  NTA31  *  NTA32  4  NTA33  •  NTAE4  points  defined  by  NTAB1 
indicated  airspeeds,  NTAB2  normal  load  factors,  NTAB3  alti¬ 
tudes,  NTAB4  weights.  These  points  are  entered  as  ordinates 
of  the  simple  graph  TABLE  which  was  defined  in  paragraph  2  of 
this  section. 

The  first  card  for  the  stress  table  contains  five  (5)  fixed 
point  rumbers  in  i5  fields  ir.  the  order 

0)  l 
f  2)  T’TAri 
(3)  NT  A:  32 
(1)  T7AE3 
(5)  NTAE4 

The  next  caru  ^)contain(s)  the  indicated  airspeeds  (floating 
point  numbers)  TABLE ( 1  >  through  TAbLE(i<TABl )  arranged  in 
Ei0.6  fields,  six  numoers  per  card. 

Tne  next  entries  are  trie  normal  loaa  factors  TABLE(NTAB1  +1, 
NU12)  (see  paragraph  2  for  definition  of  arguments)  ar¬ 
ranged  in  E10.3  fields  six  numbers  per  card. 

Next,  the  card(s)  that  contain  the  altitudes  TABLE (NN’i 2+1 ) 
through  TABLE(NN13)  arranged  in  E10.3  fields,  six  numbers 
per  curb  are  entered  in  order. 


The  next  entries  depend  on  the  number  NTS. 

If  NTB  =  1  the  card(s)  that  contain(s)  the  weights  TABLE ( NNl 3+1 ) 
through  TABLE(NN14)  arranged  in  E10.3  fields,  are  entered  with 
six  numbers  per  card. 

The  next  card(s)  contain(s)  the  stresses  TABLE(NN14+1 )  through 
TABLE(NF)  arranged  in  E10.3  fields,  six  numbers  per  card.  The 
ordering  of  the  stresses  in  this  entry  is  defined  in  paragraph 
2  of  this  section . 

If  NTB  =  2  a  card  is  entered  that  contains  the  reference  weight 
WTTB3  in  an  E10.3  field. 

The  next  card(s)  contain(s)  the  stresses  TABLE(NN1 3+1 )  through 
TABLE(NF)  arranged  in  E10.3  fields,  six  numbers  per  card. 

The  number  NF  must  not  exceed  2000.  The  ordering  of  the  stresses 
in  this  entry  is  defined  in  paragraph  2  of  this  section. 

(Note  that  NTAB4  =  1  for  this  case.) 

This  completes  the  stress  table 

If  I  =  2  then  the  entry  is  the  VGH  histogram  table  where  there 
are  NP24  »  NT421  *  NT422  '  NT423  ’  NT424  regions  defined  by  NT421 
indicated  airspeed  intervals,  NT422  normal  load  factor  intervals, 
NT423  altitude  intervals,  and  NT424  weight  intervals. 

The  first  card  for  the  VGH  histogram  table  contains  five  fixed 
point  numbers  in  15  fields  in  the  order 

(1)  2 

(2)  NT421 

(3)  NT422 

(4)  NT423 

(5)  NT424 

Following  this  card  are  the  card(s)  with  the  indicated  airspeeds 
(floating  point  numbers)  VII(l)  through  V 1 1 ( NT421 )  arranged  in 
El 0.3  fields,  six  numbers  per  card. 

The  next  card(s)  contain  the  normal  load  factors  NZ  1(1)  through 
NZI(NT422)  arranged  in  E10.3  fields,  six  numbers  per  card. 
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Next  are  the  card(s)  that  cor  tar.  tne  altitudes  Hl(l) 
th**o-igh  H!(*!T423)  arranged  in  LIU.;'  tv-  or,  .  r<  number-- 
per  curd. 

The  weight  entries  W  1(1)  through  wI(NT4?4)  at  raruvd  in 
HO. 3  field,  six  lumbers  p«r  card,  are  next. 

The  final  cerd(s‘  in  the  VGH  hisWui,  v  iecf  arc  the  loan 
occurrences  in  regions  defined  t  -e  mt'cafed  airspeed'., 
no-mal  load  factors,  altitu  es,  and  weights.  If  a  in 
[1 ,  NP]  then  tiese  entries  arc-  1  (,  V 1 1  ( 1 ) ,  NZ I  ( 1 ) ,  MIN),  '.11(1 ) ' 
through  ya(VII  (NT4?1 ) ,  NZI(\'T4?2),  !-I(';T4"3J,  W I  ( N  T  ^  24 ) 
arranged  in  E1C.3  -cie1ds,  si;  umber,  per  card.  If  i  is  in 
[1,NT4?1],  j  is  in  (.1  ,NT4?2J,  -.  is  n  L 1  ,N  142.-  j ,  and  m  is  in 
[1.NT424]  then  the  stress  that  corresponds  to  VII (i), 

N7-T(j),  LICO,  WI(ir)  is  the ( (m-1 )  ‘  '11431  '  NT422  *  NT483  + 
(k-1)  ‘  NT421  *  NT422  +  (j-1)  *  MT421 1  i ) th  entrv  on  these 
cards.  The  number  NT421  +  NT422  +  ,\!T423  +  NT424  *  NT421  • 

NT422  *  NT423  *  MT424  must  not  exceed  2000. 

6  SAMPLE  PROBLEM 

A  sample  run  is  presented  for  the  purpose  of  acquainting  the 
user  with  the  input  data  cards  and  the  outwit.  The  data  used 
does  not  represent,  any  particular  aircratt  or  usage.  It  is  as¬ 
sumed  that  two  control  points  are  sufficient  in  this  case  to  de¬ 
fine  the  full  scale  aircraft  fatigue  spectrum.  The  innut  cards 
are  as  follows: 

ICO  5  0  0  0  2  I  .3  3U  1973  23  1  2  2  21 

2  2  2  2  0  0  1 

CHECK  OUT  RUN  FOR  SPECA  PROGRAM 

V  G  H  DATA  IN  TABLE  CONTROL  POINT  NUMBER  1  (REVISION  2) 

1  6  1 


650.0  7.0 

40000.0 

36000.0 

1  .0 

500.0 

101  123  1 

20000.0  22000.0 

24000.0 

26000.0 

28000.0 

30000.0 

32000.0  34000.0 

36000.0 

38000.0 

40000.0 

42000.0 

44000.0  46000.0 

48000.0 

50000.0 

52000.0 

54000.0 

56000.0  58000.0 

60000.0 

62000.0 

64000.0 

?m  21 2  1 

0.05  0.10 

0.20 

0.25 

0.30 

0.40 

0.50  0.60 

0.70 

0.80 

0.90 

0.95 

301  302  1 

20000.0  35000.0 

r 


TW 


326  3 

25000.0 


27  1 

30000.0 

1 


300.0 

500.0 

600.0 

3.0 

6.0 

8.0 

5000.0 

20000 . 0 

35000.0 

25000.0 

30000.0 

35000.0 

20000.0 

22000.0 

25000.0 

42000.0 

43000.0 

46000,0 

55000.0 

57000.0 

58000 . 0 

17000.0 

18000.0 

20000.0 

32000.0 

34000.0 

35000.0 

51000.0 

53000.0 

54000.0 

1 5000.0 

17000.0 

18000.0 

26000.0 

28000.0 

29000.0 

4 1 000 . 0 

42000.0 

44000.0 

24000.0 

26000.0 

29000.0 

46000.0 

47000.0 

50000.0 

59000.0 

61000.0 

62000.0 

21000.0 

22000.0 

24000.0 

36000.0 

38000.0 

39000.0 

55000.0 

57000.0 

58000.0 

19000.0 

21000.0 

22000.0 

30000.0 

32000.0 

33000.0 

45000.0 

46000.0 

48000.0 

27000.0 

29000.0 

32000.0 

49000.0 

50000.0 

53000.0 

62000.0 

64000.0 

65030.0 

24000.0 

25000.0 

27000.0 

39000.0 

41000.0 

42000.0 

58000.0 

60000.0 

61000.0 

22000.0 

24000.0 

25000.0 

33000.0 

35000.0 

36000.0 

48000.0 

49000.0 

51000.0 

2 

3 

3 

3 

3 

350.0 

450.0 

550.0 

4.0 

5.0  6.0 

10000.0 

20000.0 

30000.0 

30000.0 

32000.0 

34000.0 

500.0 

1000.0 

600 . 0 

1000.0 

1100.0 

600.0 

700.0 

770.0 

7700.0 

100.0 

4000.0 

600.0 

1000.0 

100.0 

1000.0 

700.0 

2000.0 

600.0 

6000.0 

1000.0 

4000.0 

600.0 

6000.0 

1000.0 

700.0 

1000.0 

6000.0 

3000.0 

600.0 

700.0 

1000.0 

2000.0 

500.0 

4000.0 

700.0 

1000.0 

500.0 

4000.0 

6000.0 

700.0 

7000.0 

7700.0 

4000.0 

2000.0 

700.0 

500.0 

3000.0 

500 . 0 

1000.0 

600.0 

3000.0 

1000.0 

500.0 

5000.0 

700.0 

4000.0 

600.0 

3000.0 

4000.0 

400.0 

1000.0 

100.0 

500.0 

5000. 0 

5500.0 

4000.0 

500.0 

1000.0 

700.0 

500.0 

5000.0 

3000.0 

600.0 

6000.0 

3000.0 

2000.0 

4000.0 

3000.0 

1000.0 

2.000.0 

200.0 

i  nn 

n  n 

0  0 

1  1 

8  30  1973  Z3 

300.0 


2  2  2  0 
CHECK  OUT  TOR  SPEC 
V  G  H  DATA  IN  TABL 
1  3 

500.0  600.0 


0  1 
,A  PROGRAM 

X  CONTROL  POINT  NUMBER 
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1  2  2  21 


3.0 

6.0 

8.0 

5000.0 

20000.0 

35000.0 

25000.0 

30000.0 

35000.0 

10000.0 

20000.0 

24000.0 

45000.0 

50000.0 

53000.0 

25000.0 

27000.0 

30000.0 

15000.0 

17000.0 

18000.0 

32000.0 

34000.0 

36000.0 

42000.0 

43000.0 

44000.0 

21 000.0 

22000.0 

24000.0 

55000.0 

57000.0 

58000.0 

30000.0 

32000.0 

33000.0 

27000.0 

29000.0 

32000.0 

62000.0 

64000.0 

65000.0 

39000.0 

41000.0 

42000.0 

22000.0 

24000.0 

25000.0 

48000.0 

49000.0 

51000.0 

36000.0 

39000.0 

40000.0 

17000.0 

18000.0 

20000.0 

51000.0 

53000.0 

54000.0 

26000.0 

28000.0 

29000.0 

25000.0 

26000.0 

27000.0 

55000.0 

56000.0 

57000.0 

36000.0 

38000.0 

39000.0 

19000.0 

21000.0 

22000.0 

45000.0 

46C00.0 

48000.0 

49000.0 

50000.0 

53000.0 

24000.0 

25000.0 

27000.0 

58000.0 

60000.0 

61000.0 

33000.0 

35000.0 

36000.0 

Based  on  this  input  the  following  output  listing  was  obtained 
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SECTION  IV 

EXAMPLE  PROBLEM  -  F-4  STRESS  SPECTRUM 
FOR  POSITIVE  LOAD  FACTORS 

The  data  base  for  this  problem  is  four  quarters  of  VGH  data 
starting  with  the  second  quarter  of  1972  and  finishing  with  the 
first  quarter  of  1973.  The  VGH  histogram  intervals  for  these 
data  are  the  following: 

Indicated  airspeed  (knots) 

150,  200,  250,  300,  350,  400,  450,  500,  550,  625,  and 
700 

Normal  load  factor  (equivalent) 

1.4,  1.8,  2.2,  2.6,  3.0,  3.8,  4.6,  5.4,  6.6,  7.8,  and  9.0 
Altitude  (feet) 

0.  1000,  2000,  5000,  10,000,  15,000,  20,000,  30,000, 

40,000,  and  50,000 

Weight  (pounds) 

37,500  (reference  weight) 

The  stress  table  was  set  up  with  the  following  indicated  airspeed, 
normal  load  factor,  altitude,  and  weight  combinations: 

Indicated  airspeed  (knots) 

175,  225,  275,  325,  375,  425,  475,  525,  575,  and  625 
Normal  load  factor 

2.4,  2.8,  3.4,  4.2,  5.0,  6.0,  7.2,  and  8.9 
Altitude  (feet) 

500,  1500,  3500,  7500,  12500,  17500,  25000,  and  35000 

Weight  (pounds) 

37,500 

The  VGH  histogram  table  was  made  up  using  all  of  the  available 
data  for  the  air-to-air  and  air-to-ground  operations  for  the  four 
quarters  without  distinguishing  the  various  F-4  models  except  that 
only  the  unslatted  configurations  were  considered.  The  numbers  of 
hours  of  data  in  each  category  and  their  corresponding  numbers  of 
positive  and  negative  load  occurrences  are  shown  in  Table  1. 
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The  number  of  stress  exceedances  per  1000  hours  for  load 
reference  station  (LRS )  180,  defined  in  Figure  1  is  shown  in 
Figure  2  through  Figure  7.  Figure  2  through  Figure  5  shows  the 
variation  from  quarter  to  quarter  of  the  VGH  data.  The  stress 
exceedance  graphs  appear  to  show  a  small  degree  of  scatter 
except  for  the  SEA  air-to-air  first  quarter  where  there  was 
an  overt  change  in  the  mission  although  it  was  still  categorized 
as  air-to-air.  The  four  quarters  of  data  are  combined  in  Figures 
6  and  7  to  show  the  differences  between  the  CONUS  and  SEA  in 
the  air-to-air  operation  and  the  air-to-ground  operation. 
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SECTION  V 


CONCLUSIONS 


The  procedure  described  in  this  report  can  eliminate  much 
of  the  uncertainty  that  can  occur  in  the  derivation  of  the  maneuver 
load  stress  spectrum.  For  new  aircraft  an  estimate  must  be  made 
of  the  VGH  histogram  to  obtain  the  spectrum.  This  estimate  can 
be  updated  during  Task  V  of  ASIP  to  derive  a  better  estimate  for 
the  operational  life  of  the  fleet.  This  procedure  can  be 
immediately  applied  to  fleet  tracking  by  computing  the  conditional 
probability  of  exceeding  a  stress  level  given  the  normal  load 
factor. 

The  application  to  full  scale  aircraft  testing  makes  use  of 
the  assumption  that  the  stress  is  matched  at  a  specified  number  of 
control  points  by  a  linear  combination  of  the  same  number  of 
balanced  loading  conditions.  This  technique  is  believed  to  be  more 
accurate  than  the  usual  process  of  a  damage  match  at  the  specified 
control  points  in  that  the  troublesome  damage  calculation  is 
eliminated.  The  stress  spectra  at  points  other  than  control  points 
are  presumed  to  be  matched  satisfactorily  by  using  representative 
loading  conditions.  It  is,  of  course,  theoretically  better  to  use 
all  points  of  the  sky  that  occur  in  the  VGH  histogram.  This, 
however,  may  be  impractical  due  to  test  equipment  limitations. 

The  procedure  as  applied  to  the  F-4  fleet  indicates  that 
in  general  the  stress  spectra  do  not  show  significant  changes 
from  quarter  to  quarter.  Also,  when  an  operational  change  is  made 
the  method  will  reflect  that  change.  When  CONUS  and  SEA  data  are 
compared  there  appears  to  be  a  reasonably  good  correlation  between 
the  spectra  generated  in  training  and  the  spectra  generated  in 
combat . 
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TABLE  1 .  F-4  VGH  DATA  SUMMARY 


PERIOD 

TYPE 

HOURS 

+  COUNTS 

-  COUNTS 

HHJQV 

2Q  72 

CONUS 

AA 

196.87 

13981 

4541 

18522 

2Q  72 

SEA 

AA 

65.84 

4C47 

562 

4709 

2Q  72 

CONUS 

AG 

251 .70 

307  23 

9432 

40155 

2Q  72 

SEA 

AG 

1248.94 

93027 

13258 

106285 

3Q  72 

CONUS 

AA 

290.22 

17470 

5259 

22729 

3Q  72 

SEA 

AA 

393.26 

30842 

6366 

37208 

3Q  72 

CONUS 

AG 

469.64 

38968 

9934 

48902 

3Q  72 

SEA 

AG 

802.74 

6644  6 

12485 

78931 

4Q  72 

CONUS 

AA 

184.92 

9404 

3138 

12542 

4Q  72 

SEA 

AA 

164.35 

7862 

1933 

9795 

4Q  72 

CONUS 

AG 

89.16 

5959 

1265 

7224 

4Q  72 

SEA 

AA 

502.09 

28254 

4819 

33073 

IQ  73 

CONUS 

AA 

123.18 

7983 

2231 

10214 

IQ  73 

SEA 

AA 

133.96 

6100 

1358 

7458 

IQ  73 

CONUS 

AG 

194.38 

17828 

4001 

21829 

IQ  73 

SEA 

AG 

933.87 

42342 

. -  , 

8229 

50571 

52 


3Q  72 
4Q  72 
IQ  73 


Exceedances  per  1000  hours 


Exceedances  per  1000  Hours 


Stress  Exceedances  per  1000  hours 


Stress  -  kips 

Figure  6.  F-4  Spectra  -  Air-to-Ai 

for  One  Year  of  VGH  Oat 


stress  -  kips 

F-4  Spectra  -  Air-to-Ground  (All  Models) 
for  One  Year  of  VGH  Data 


APPENDIX  -  SPtCA  PROGRAM  LISTING 

The  listing  given  below  is  a  FORTRAN  extended  language 
routine.  This  listing  contains  all  of  the  statements  for 
the  version  described  in  the  Introduction  of  this  report. 
Section  3.3  gives  a  brief  description  of  each  of  the  sub¬ 
routines  in  this  listing. 
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PROGRAM  FOR  COHPUlIKG  AIRCRAFT  INTERNAL  LOAO  PRCBABILITt 
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